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ABSTRACT 

We investigate the spectral variability of the Seyfert galaxy Fairall 9 using almost 6 years of monitoring 
with the Rossi X-ray Timing Explorer (RXTE) with an approximate time resolution of 4 days. We 
discover the existence of pronounced and sharp dips in the X-ray flux, with a rapid decline of the 2-20 
keV flux of a factor 2 or more followed by a recovery to pre-dip fluxes after ^ 10 days. These dips skew 
the flux distribution away from the commonly observed log-normal distribution. Dips may result from 
the eclipse of the central X-ray source by broad line region (BLR) clouds, as has recently been found in 
NGC 1365 and Mrk 766. Unlike these other examples, however, the clouds in Fairall 9 would need to 
be Compton-thick, and the non-dip state is remarkably free of any absorption features. A particularly 
intriguing alternative is that the accretion disk is undergoing the same cycle of disruption/ejection as 
seen in the accretion disks of broad line radio galaxies (BLRGs) such as 3C120 but, for some reason, 
fails to create a rclativistic jet. This suggests that a detailed comparison of Fairall 9 and 3C120 with 
future high-quality data may hold the key to understanding the formation of rclativistic jets in AGN. 
Subject headings: galaxies: individual(Fairall 9) - X-rays: galaxies - galaxies: nuclei - galaxies: Seyfert 



1. AGN VARIABILITY 

X-rays are emitted from the inner accretion disks 
around both stellar-mass and supermassive black holes. 
This gives us our most effective probe to date of accre- 
tion physics close to black holes as well as the immediate 
environments of black holes. Studies of X-ray spectral 
variability are a particularly powerful diagnostic. For ex- 
ample, combined spectral and timing analyses have led 
to the dis covery of distinct hardness /flux states in X-ray 
binaries (jTanaka fc Shibazakil 119960 which almost cer- 
tainly correspond to distinct accretion modes. Spectral 
variability also facilitates the decomposition of a com- 
plex spectrum into physical components that would be 
degenerate in a single spectrum. 

Detailed studies of long term (year-f ) X-ray variabil- 
ity in active galactic nuclei (AGN) have been possible 
since the launch of the Rossi X-ray Timing Explorer 
(RXTE) but a relatively small n umber of sources have 
been targeted for detailed study |Markowitz et al.ll2003l : 
iRivers. Markowitz fc Rothschild 1201 It ). These studies 
are crucially important as they probe timescales com- 
parable to the viscous timescale on which the true mass 
accretion rate would be expected to vary. More typical 
X-ray observations of active galactic nuclei (AGN) probe 
timescales of hours-to-days, characteristic of dynamical 
or thermal timescales but much shorter than the viscous 
timescales. Long-term studies are also important for un- 
covering rare/transient phenomena that may give unique 
windows into the physics and structure of AGN. 

In this Letter, we report the discovery of an un- 
usual pattern of X-ray variability in the nearby Seyfert- 1 
galaxy Fairall 9 (z = 0.047) using the long-term mon- 
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itoring data from pointed RXTE observations. Rather 
than displaying the usual "flary" lightcurvc with a 
log-normal flux distribution jGaskcU & Klimck 200^ 
lUttlev. McHardv fc Vaughanl 1 2005( 1 . Fairall 9 shows 
long-timescale modulations punctuated by strong and 
sudden dips. This behavior was apparent but not com- 
ment ed upon when the m onitoring was first presented 
by M arkowitz et al.l ()2003( ) . Using the latest background 
models and calibration files, we analyze all available 
RXTE data and confirm the presence of these dips in 
the RXTE lightcurve. We discuss their possible origins, 
including the possibility that they correspond to eclipses 
of the X-ray source by clouds in the broad line region 
(BLR). Alternatively, the dips may result from a gen- 
uine accretion disk instability, suggesting an interesting 
parallel between Fairall 9 and broad line radio galaxies. 

2. OBSERVATIONS 

The observations considered in this analysis are all 
those of Fairall 9 with pubhc RXTE /FCA data, a to- 
tal of 744 pointings spanning the period from late 1996 
until early 2003. This gives an average time resolu- 
tion of about 4 days with an average pointing length of 
1.1 ks. Only data from the top layer of proportional 
counter uni t 2 were use d . Th e data were reduced as 
outlined in iWilms et al.l (|2006f ). using Heasoft 6.10. 
For each pointing, we produce a background-subtracted 
PGA spectrum using the appropriate epoch-dependent 
background/response files. The count rates quoted be- 
low were obtained from the spectra for each pointing. 
We do not consider any intra-pointing variability in this 
work. 

3. RESULTS 

3.1. Light curves and confirmation of dipping 

The long-term light curve in the 2-10 keV and 10- 
20 keV bands is plotted in Figure[T](left). This light curve 
displays the unusual behavior already noted in the intro- 



Lohfink et al. 



JD-2450000 

400 700 1000 1300 1600 1900 2200 2500 2800 



JD-2450000 

1800 



-JTTTTJTTTTjn 

lOkeV 



'■^* ^^I&saM^'^'P 



y '1 , '111 J ,« " /' 



1.00" 




1.00 



S^ 0.10 



0.01 



2-lOkeV 



I , I 



I V 



1 0-20 keV 



'I, III 



III,, I |l 



'l l,l' , I I,l|t 



|l I I I I II III 'il , 'ii 'M'. 

' I'lll III, I I I'l I I I I V I 



1997 1998 1999 2000 2001 2002 



Jun Jul Aug Sep Oct Nov Dec Jan Feb Mar 
200(5 2001 



E 40 



l_l 

, ^ 


s 


u 


JL 


vJ 



0.1 1.0 

net rate [2-10 keV] 
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Figure 1. Left : Long term RXTE-PCA light curves for 2-10 IceV (top panel) and 10— 20keV (lower panel). Right : Zoom-in on the 
RXTE-PCA light curves for mid 2000 until 2001 April. 

To further explore the nature of the variability, we per- 
form flux resolved spectroscopy. We choose to use the 
2-5 keV flux as our flux discriminator since this will be 
dominated by the primary power-law rather than a soft 
excess component or reflection. We define a "high state" 
as F2_5kcV > 1 X 10"^-'^ ergss"^ cm~^ and, by summing 
the data for all pointing which exceed this flux, form a 
"high state spectrum" . Similarly, we define a "low state 
spectrum" (F2-5keV — (0.5-1.0) x 10~^^ ergs s~^ cm~^). 
A flux of 10^^^ ergs s""'^ in 2-5 keV corresponds to a 2- 
lOkeV count rate of about 1.45 counts per second. A 
"dip state spectrum" , with a dip defined as a count rate 
drop of a factor of 1.5 in 2-10 keV, with respect to the 
rate right before the dip and the two pointings following 
the drop, is also defined. The pointings considered to be 
part of a dip were excluded from high and low states. 

Previous studies have shown that, to first order, the 
spectrum can be well described by a continuum, orig- 
inating from a comptonizing corona, and cold reflec- 
tion. We model th is scenario using the pexmon model 
(jNandra et al.|[2007l ). which consists of a power law con- 
tinuum plus the cold reflection continuum and associ- 
ated iron Ka/K/3 fluorescent lines (with a self-consistent 
strength). The quality of the data does not allow con- 
straints on all parameters in the model, therefore the 
abundances are assumed to be solar and the inclination 
of the reflector is fixed to 60 degrees (the most prob- 
able value, assuming a random orientation). We also 
include cold Galactic absorption with a column density 
iVn ~ 3.1 X lO^'^cm "^, modeled b y TBneiQ a newer ver- 
sion of TBabs (Wilms. Allen fc McCrav.2000.) with cross 
sections set to vern and abundances set to wilm. 

Table[T]shows the spectral fitting results for these three 
flux-sorted spectra. The spectra are well described by 
the model, see Fig. [31 except for the high state spec- 
trum where a residual line in the iron K region remains. 
This feature can be modeled with an additional iron line 
at the redshift of Fair all 9, with an equivalent width of 
35 ± 8 eV. The uncertainties in the energy of this addi- 
tional line span both neutral and H-like iron. Beyond 



duction; long timescale flux modulations punctuated by 
short and intense dips. These dips consist of a rapid flux 
decline by a factor of 2-4 in 5-10 days and then a re- 
covery to the pre-dip level with an entire dip-duration of 
10-20 days (Fig. [T] right). The existence of these dips in- 
troduces a low-flux tail to the flux distribution, skewing 
it strongly away from log- normal (Fig. [2]). To the best of 
our knowledge, these unusual dips in Fair all 9 have not 
been commented upon previously. 

On longer timescales the source is experiencing a flux 
decline from early 1999 reaching its low in the fall of 
2000 and then recovering to a higher flux value towards 
the middle of 2001. There is clearly a greater propensity 
to dip during the lower-flux state. The average flux in 2- 
lOkeV for the object is (2.16±0.63) -lO"" ergss'^ cm'^. 
Both the long term and short term variability seem to be 
of similar nature in the two different bands. In particu- 
lar, while the dips are most clearly seen in the 2-10 keV 
light curve due to the higher signal-to-noise, they are also 
present in the 10-20 keV light curve. 



3.2. Spectral Evolution 



' http://pulsar.sternwarte.uni-erlangen.de/wilms/research/tbabs/ 



X-ray dips in Fairall 9 
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Figure 3. Spectra (high state [dots], low state [crosses], dip state 
[stars]), model [solid line] and residuals for the PCA summed flux 
resolved spectra for the fit with the eold refieetion model. 



Table 1 

Spectral Parameters for eold refieetion fits, see text for details on 
the model set-up. 
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0.5«-J 

1 1+0.3 
^•^-0.2 

2 1+1-3 


2.7-10"" 
2.7- 10"" 
1.8-10"" 
1.2- 10"" 


73.3/42 
16.5/41 
38.9/42 
32.0/42 



that required by the reflection component, an additional 
hne -with this strength is rejected by both the lo-w-state 
and dip-state spectra. Thus, it appears that the addi- 
tional iron line component is exclusively displayed by the 
high state spectrum. We comment on the nature of this 
spectral feature in Section 4. 

In agre ement -with previou s studies of Seyfert galaxies 
(e.g., see lChiang et al.l 120000 . the continuum softens (F 
increases) from the lo-w- to the high-state (Fig. |3]). We 
also see a reflection fraction -which is inversely propor- 
tional to the continuum flux; i.e., the X-ray reflection is 
consistent -with being a constant flux contribution. The 
most obvious interpretation is that a signiflcant compo- 
nent of the reflection is from a distant structure. Given 
that the modulation bef-wen high- and low-flux states oc- 
curs on timescales of ~ 1 year, we infer that the distant 
reflector is situated at least a light year from the cen- 
tral X-ray source. We note that X-ray reflection from 
the inner accretion disk is also expected to be pr esent 
(jSchmoU et al.ll2009l: lEmmanoulopoulos et al.ll2QllD and, 
given the error bars on R, its presence is consistent with 
these data. 

Although flux resolved spectroscopy is a po-werful tool, 
it does not offer the possibility to study the actual spec- 
tral evolution in time. Higher time resolution can be 
achieved by looking at X-ray flux-flux plots in different 
bands. Figure |3] sho-ws the 2-5keV/5-10keV flux-flux 
plot; a positive correlation is apparent, although it is 
clearly non-linear as can be seen by a comparison -with 
a spectral variability model possessing a flux variable 
po-wer la-w with a constant hard component (i.e. constant 
cold reflection). At high fluxes, the non-linearity appears 
dominated by a genuine softening of the primary contin- 
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Figure 4. Flux-flux plot for the 2-5 keV and the 5-10 keV band 
binned for plotting, showing the non-linearity of the flux-flux rela- 
tion, as well as model lines for a constant hard component with a 
fiux- variable power law for a set of photon indexes. 

uum. This was already seen in earlier studies for Fairall 9 
and other Seyfert 1 galaxies ("IVIarkowi tz fc Edelsonll2001l : 
[Sobolewska & Papadakis 2009). At the lowest fluxes 
(corresponding to the dip state) there is a change in 
the nature of this non-linearity, with the lowest few 2- 
5keV flux points having an almost constant 5-10 keV 
flux. This signals a significant spectral hardening during 
the deepest parts of the (2-5 keV) dips. However, the 
spectrum is never dominated by reflection. 

We do note that this non-linearity calls into question 
the use of y-intercepts in extrapolated flux-flux diagrams 
to determine constan t components in AGN spectra (e.g., 
see lNodaet~aI|[20Tl . 

4. DISCUSSION AND CONCLUSION 

In this letter, we studied the X-ray variability of 
Fairall 9 on timescales of a few days to five years. On 
the long timescales studied by the flux resolved spec- 
troscopy, a significant fraction of the refiected X-rays ap- 
pear to maintain a constant normalization, suggesting 
that potentially a significant fraction of the X-ray re- 
flection might originate from a structure at least a light 
year across. The excess iron emission line (above and 
beyond that associated with the reflected component) 
in the high-state spectrum suggests that, in this state, 
there might be an appreciable increase in the amount of 
Compton-thin material intercepting the primary X-ray 
continuum. This material must be out of our line of 
sight since we see no iron absorption edge in the high- 
state spectrum. A possible identiflcation for these feature 
is iron line emission (fluorescence or radiative recombi- 
nation) from a radiatively-driven disk wind that forms in 
the high-flux state. The spectral resolution of the PCA 
prevents us from constraining the charge state of the ad- 
ditional iron line and hence the ionization state of this 
additional material. 

The most unusual flnding is the discovery of strong and 
sudden dips in the observed X-ray flux from Fairall 9. In 
the rest of this section, we discuss the nature of these 
dips. There are two fundamentally different interpreta- 
tions; transient absorption events (i.e. eclipses) of the 
inner accretion disk, or a genuine turning-off of X-ray 
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emission from the inner disk. 

We begin by addressing the absorption possibihty. At 
least upon the immediate onset of a dip, both the 2- 
lOkeV and 10-20 keV fluxes drop by similar magnitudes 
requiring that the absorber be Compton-thick {N^ > 
few X 10^"* cm~^). Suppose that the radius of an absorb- 
ing "cloud" has radius r = xrg (where rg = GM /c^ is 
the gravitational radius) and it is situated R = Xvg from 
the black hole. Assuming Keplerian motion, the veloc- 
ity of the cloud is v — X~^/'^c and hence the eclipse 
of the inner disk by the cloud with take a duration 
At ^ 2r/v — 2xX^/' ^{rg/c). For a black hole mass 
of M K. 2> X IQ^Mq (jPetersonet al.ll2004| ). this gives 
At - 3000a;Xi/2s. The observed 2-10 keV dips have 
duration Ai « 2 x 10^ s, implying that xX^/"^ w 700. 
Assuming that the X-ray source has a size of rx = lOr^, 
we must have a; ^ 10 (r > 5 x 10^^ cm) in order for 
the cloud to be able to block the source and create a 
strong dip, implying X < 5000 and v > 4300 km s~^. 
Thus, the absorbing clouds are consistent with being at 



BLR distances. For a cloud with radius 



10 



15, 



to be Compton-thick, it must have a density exceeding 
n ~ few x lO^cm"'^. Thus, these absorbing clouds also 
hav e a density consistent w ith that of BLR clouds (e.g., 
see iKorista fc G"oadll2000l ) 

X-ray eclipses by putative BLR clouds have been 
found recently in som e other AGN, notably NGC 136 5 
(jMaiolino et al.ll2010D and Mrk 766 (Risa liti et alll2011[ ). 
There are some important differences between the 
Fairall 9 events and the events in these other sources. 
Firstly, in NGC 1365, even the non-eclipse spectrum 
shows significant absorption. By contrast, the non-dip 
state of Fairall 9 is notable for being so clean, having 
no discernible cold or ionized absorption seen in any 
of its previous ^gC^ A (Rcvnolds 1997), XMM-New ton 
(iGondoin et al.l 120011: Emnianoulopoulos et~al]l2011[ ) or 



Suzaku (jSchmoll et al.ll2009l ) observations. Secondlv. the 
eclipsing clouds in NGC 1365 and Mrk 766 are gener- 
ally Compton-thin with strong evidence that the col- 
umn density spikes up and then steadily decreases dur- 
ing the event. This l ed to the notion of BLR "comets" 
(jMaiolino et al.l 120101 ) with high column density heads 
followed by lower column tails. By contrast, the Fairall 9 
clouds must be Compton-thick, at least at the beginning 
of the dip event. However, a "comet" scenario would pre- 
dict a Compton-thick to Compton-thin transition, lead- 
ing to a recovery of the hard band light curve before the 
soft band dip ends. The signal to noise of the current 
data is insufficient to detect such an early hard-band re- 
covery. Future studies with Suzaku. NuSTAR or Astro-H 
are required to test these predictions. 

The presence of well-defined Compton-thick clouds in 
an otherwise extremely clean environment raises ques- 
tions about cloud confine ment and, more general l y, the 
nature of thes e clouds (jKallman fc Mushotzkvl 119851 : 
iSnedden fc Gas kcll 200'3). A detailed discussion of the 
physics of such clouds is beyond the scope of this Letter. 
We do however note that the properties of these clouds 
(Compton-thick and r ^ 10^^ cm) are consistent with the 
irradia ted envelopes of circumnuclea r stars (i.e., "bloated 
stars" ■ lEdwa7d^ll980l : lPenstonlll98 8). Bloated stars have 
been ruled out as a viable model for the broad emis- 



sion lines the mselves on the basis of the smooth ness of 
the Ha profile (jArav et al.lll997HLaor et al.ll2006D . How- 
ever, especially in high-mass systems such as Fairall 9, 
bloated stars may nevertheless be present and could read- 
ily eclipse the X-ray source producing dips. 

Could the dips be intrinsic to the accretion process, 
i.e., signal a genuine "turning-ofF' of X-ray emission from 
the inner accretion disk? It will be useful to delineate a 
few characteristic timescales for the inner regions of this 
AGN. The light crossing time of the inner X-ray emit- 
ting regions of the accretion disk (say, out to r = lO^g) 
is iic ~ 30 ks, and the dynamical timescale of the disk at 
r = lOvg is tdyn ~ 50 ks. Assuming an angular momen- 
tum transport parameter a ~ 0.1, the thermal timescale 
at this radius is Uh = tdyn/ct w 500 ks. Finally, we con- 
sider the viscous timescale, ivisc = Uh/ih/r)'^ where h 
is the geometric thickness of the accretion disk. With 
an average 2-10 keV fiux of 2.2 x 10~^^erg s/sec/cm^ and 
using a bolometric correction of 50 from iMarconi et al.l 
(|2004D . we find an Eddington ratio of £ w 0.15 for 
Fairall 9. For such an Eddington ratio (and assuming 
a radia tive efficiency of ?? = 0.1), s tandard accretion disk 
theory (jShakura fc Sunyaevlll973l ). the inner disk is radi- 
ation pressure dominated and has a geometric thickness 



{SC/2r^)rg 



2.3 To- Thus, at r 



10 Tg, we have 



h/r = 0.23 and a viscous timescale of tvisc ~ 9.4 Ms (i.e. 
approximately 110 days). 

The most obvious (albeit dramatic) scenario is one 
in which the innermost parts of the radiatively-efficient 
optically-thick accretion disk are destroyed/ejected by a 
dynamical or thermal instability. This is hypothesized 
to occur in radio-loud broad line radio g alaxies (BLRGs) 
as part of the process of jet formation (iMarscher et al.l 
2002; Chattcri ee et all 120091 : King et al.ll20lH ). In the 
BLRGs, a distinct dip in the X-ray flux precedes large 
radio flares and the creation of a new superluminal knot. 
Since, as with Seyferts, the X-ray emission in BLRGs 
is als o dominated by the corona of the inner accretion 
disk (Mar shall et al.l [20091) . this X-ray/radio connection 
conclusively demonstrates a link between changes in the 
inner accretion disk structure and powerful jet ejection 
events. In fact, it is instructive to compare directly 
Fairall 9 with the BLRG 3C120. 3C120 has a slightly 
lower-mass black hole, 6 x 10^ Mc^, for 3C120 c ompared to 
3 X 10^ A/0 for Fairall 9 (jPeterson et al.ll2004[ ). However, 
they possess a very similar Eddington ratio, Cpg « 0.15 
compared to £30120 ~ 0.11, derived using the bol o metri c 
luminosity for 3C120 from iVasudevan fc FabianI (|2009| ). 
Thus, their accretion disks might be expected to be in 
rather similar regimes of behavior. It is possible that 
the X-ray dips seen in Fairall 9 have the same physical 
origin as the ejection-related X-ray dips seen in 3C120 
but, for some reason, these disk disruption events lead 
to the creation of powerful jet outflows in 3C120 but not 
in Fairall 9 (indeed, Fairall 9 has yet to be detected in the 
rad io band, with a n upper limit of < 36mJy at 1.4 GHz, 
see iWhittld 11993 ). Either the magnitude/direction of 
the black hole spin or the magnetic field structure in 
the inner disk may be the determining factor in decid- 
ing whether disk disruptions generate relativistic out- 
flows. For Fairall 9, the black hole spin ha s already been 
meas ured and was found to be moderate (|Schmoll et al.l 
l2009f ). Observations of Galactic black hole X-ray bina- 
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ries support this idea of viewing the dips in the context 
of jet ejections, as some were observed to possess very 
faint jets, undete ctable in case of an AGN ([Gaho et aLl 
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